First direct observation of the Van Hove singularity in the tunneling spectra of cuprates 
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In two-dimensional lattices the electronic levels are unevenly spaced, and the density of states (DOS) 
displays a logarithmic divergence known as the Van Hove singularity (VHS). This is the case in particular for 
the layered cuprate superconductors. The scanning tunneling microscope (STM) probes the DOS, and is 
therefore the ideal tool to observe the VHS. No STM study of cuprate superconductors has reported such an 
observation so far giving rise to a debate about the possibility of observing directly the normal state DOS 
in the tunneling spectra. In this study we show for the first time that the VHS is unambiguously observed 
in STM measurements performed on the cuprate Bi2Sr2Cu05+5 (Bi-2201). Beside closing the debate, our 
analysis proves the presence of the pseudogap in the overdoped side of the phase diagram of Bi-2201 and 
discredits the scenario of the pseudogap phase crossing the superconducting dome. 
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INTRODUCTION 

The scanning tunneling spectroscopy (STS) is known to 
probe the density of states, and a logarithmic peak is 
therefore expected in the tunneling spectrum of materials 
having a VHS near the Fermi level. How^ever, there is a 
strong theoretical debate around the possibility to observe 
such features in tunneling spectra. The compelling absence 
of such a peak in the tunneling spectra of the cuprates has 
lead to claims that the normal-state DOS w^ould in fact be 
canceled like just as in ideal planar tunnel junctions. ^'"^ An 
alternative is that the VHS peak is suppressed due to the 
interaction of electrons with collective excitations.^"^ 

Among the cuprates, Bi-2201 has a combination of two 
properties required for an optimal STS investigation of the 
VHS: it is nearly two-dimensional, and has a low maximal 

of around 12 K giving access to the normal state at low 
temperature. Producing large and pure non cation-doped 
single crystals of Bi-2201 is difficult and consequently, this 
compound has attracted little attention. We succeeded in 
growing such crystals and tuning their macroscopic doping 
level,^ allowing doping dependent STS studies. 

In this study, we demonstrate that the VHS is directly seen 
in the tunneling spectra of Bi-2201. We observe that all the 
spectra can be fitted by a model involving a VHS and a d- 
wave BCS gap. This suggests that the VHS shows up in the 
low-T^ Bi-2201 due to a weak coupling to collective modes. 
Moreover, our measurement at a temperature above of 
a heavily overdoped Bi-2201 sample clearly exhibit a gap 
feature that we interpret as the signature of the pseudogap 
phase. 



RESULTS 

STM data for overdoped Bi-2201 samples. Low- 
temperature spatially-resolved STS of an as-grown sample 
with = 11 K and a transition width AT^ = 4.5 K are 
depicted in Fig. 1. Most of the surface presents spectra 
with a clear superconducting gap in the 15-45 meV range. 
However, extended regions covering 11% of the sample ex- 



hibit qualitatively different spectra, without a gap but with 
a single sharp peak 10-40 meV below the Fermi energy. A 
broadly peaked background has been seen earlier in Bi-2201 
but only at high temperature.^ The observation of this new 
type of spectra at low temperature, which we attribute to 
the VHS, is our main finding. To each spectrum we assign a 
characteristic energy which is either the peak to peak gap 
value (if two peaks are detected) or the energy position of 
the single peak. The transition between the two characteris- 
tic kinds of spectra is remarkably gradual and reproducible, 
both spatially (Fig. Id) and when the spectra are sorted 
according to the value of E^ (Fig. la). The general trend is 
that as the height of the peak decreases, a coherence peak 
emerges at positive energy signaling the opening of the gap, 
and the negative-bias peak shifts to larger negative energy. 
The nanometer-sized red to yellow areas corresponding in 
Fig. lb to single-peak spectra are all surrounded by green 
regions which have small gaps. 

We have then acquired E^-maps and the corresponding 
histograms for two overdoped samples (ODIOK and 0D7K; 
Fig. 2) . Comparing the histograms reveal a steady increase 
in the number of single-peak spectra with doping, from 
11% at r, = 11 K to 28% at = 10 K and 48% at T, = 
7 K. The mean peak position also shifts towards the Fermi 
energy with values —28, —26, and —21 meV at T^, = 11, 10, 
and 7 K, respectively, as expected in a rigid-band picture. 
Along with the increased proportion of single-peak spectra 
with increasing overdoping, we also observe a decreased 
proportion of wide-gap spectra that are excluded from the 
histograms, as one or both coherence peaks are missing. 
These spectra (depicted as black pixels in the E^-maps) 
have a proportion of 28% in AGllK, 6% in ODIOK, and 
5% in 0D7K. Much of this decrease is an effect of the post- 
annealing treatments which favor a more homogeneous 
distribution of oxygen, as confirmed by the drop of the 
superconducting transition width AT^. 

Finally we have studied heavily overdoped sample in 
which has been reduced to less than 1.8 K. The propor- 
tion of single-peak spectra increases to ~ 65% and some 
of these spectra even exhibit a peak at energy above Ep 
(Fig. 3a). This suggests that a transition from a (tt, tt)- 
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FIG. 1. Low-temperature (5 K) STS spectroscopy of as-grown 
Bi-2201 (AGllK). (a) Three dimensional overview of the data: 
each curve is the average of 60 spectra sharing the same Ep value 
(Ep is the peak position for single-peak spectra, and half the peak- 
to-peak distance for two-peak spectra); the curves are sorted by 
increasing Ep. Note the smooth transition between completely 
different curves at extremal values of Ep. (b) E^-map of a 14 x 
14 nm^ area. The white bar corresponds to 2 nm. The red and 
yellow regions correspond to negative Ep (single-peak spectra) 
while in green, blue, and purple regions the spectra exhibit a gap 
of width Ap = Ep. Red-yellow spots are surrounded by green 
regions where the smallest gap are measured (~ 20 meV) . Black 
patches are composed of spectra (~ 4000 in total) in which no 
gap could be extracted due to a missing left peak (or sometimes 
no peak at all); these spectra are not included in the histogram 
of panel c. (c) Distribution of Ep values. The dashed lines are 
Gaussian fits. The average peak position and average gap values 
are —28 and +29 mey respectively, (d) Spatial evolution of the 
spectra on the path depicted by the arrow in b. The path goes from 
a red region (top-most single-peak spectrum) to a black region 
where the left coherence peak is barely visible. Notice the strong 
similarity with the global evolution shown in a. 



centered hole-like to a (0, 0)-centered electron-like Fermi- 
surface takes place in the overdoped region near the border 
of the superconducting dome. Fits to ARPES data in cation- 
doped Bi-2201 have led to the same conclusion. As seen 
in Fig. 3b, as the gap vanishes and the peak moves to posi- 
tive energy, the spectra develop a remarkable electron-hole 
asymmetry with a drastic suppression of spectral weight 
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FIG. 2. Evolution of the Bi-2201 spectroscopic properties with 
overdoping. The top-left panels show the f^-maps (Ep is the 
peak position for single-peak spectra, and half the peak-to-peak 
distance for two-peak spectra) for a, an overdoped sample with 
= 10 K and AT, = 2.5 K (ODIOK), and for b, an overdoped 
sample with = 7 K and AT, = 1.5 K (0D7K). The white bars 
correspond to 5 nm. The distributions of Ep values presented in the 
histograms were normalized and compared with the histograms 
obtained in the least overdoped AGllK sample (dashed lines). 
The top-right panels are simultaneously acquired topographies: 
white regions are physically higher than black ones. They show 
the superstructure common to all Bi-based cuprates.^^'^^ Especially 
in ODIOK, the smallest gaps (green) are preferentially located on 
the minima of the superstructure. A similar correlation between 
gap size and the superstructure has been reported in Bi-2212^^ 
and Bi-2223.^^ Our measurements further show that the preferred 
location of single-peak spectra (red-yellow) is also the minimum 
of the superstructure, suggesting that the local doping level itself 
is modulated. 



for occupied states. This trend mirrors the common phe- 
nomenology observed in all cuprates with increasing gap, 
namely an excess of spectral weight at negative energy. The 
data of Fig. 3 were collected well above the bulk of the 
sample. It is therefore natural to interpret the series of spec- 
tra in Fig. 3b as a transition from the non-superconducting 
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FIG. 3. Spectral properties of strongly overdoped Bi-2201 
(Tc < 1.8 K) measured at 5 K. (a) Distribution of Ep values 
for single-peak spectra (red part of the histogram corresponds to 
the position of the single peak) and gapped spectra (green part of 
the histogram), (b) Representative spectra showing single peak at 
positive energy (top), single peak at negative energy and gapped 
structures. This series was obtained by ordering the spectra ac- 
cording to the position of their maximum, (c) Interpretation of b 
as a transition from the pseudogap to the normal state. T*, the 
frontier between normal and pseudogap states. 




Theoretical modeling of the tunneling spectra. In the 

first STS studies of Bi-2201, the shape of the normal-state 
spectra was not addressed. ^'^^ More recently sharp features 
in the tunneling characteristics have been reported in lead- 
doped Bi-2201. The observed structure was a narrow peak 
very close to the Fermi level (±2 meV) superimposed on 
the superconducting gap. Such spectra were found at iso- 
lated locations on the surface and, based on earlier works 
in Bi-2212,^'^'^^ they were convincingly explained as na- 
tive impurity resonances. The single-peak spectra in our 
measurements show no trace of a superimposed supercon- 
ducting gap, cover extended regions of the surface, and have 
their maxima at a doping-dependent energy up to 40 meV 
below the Fermi level. In order to support our interpreta- 
tion that the single-peak spectra are direct manifestations 
of the VHS in the DOS, we have modeled the dl/dV curves 
obtained in our samples. Following previous analysis for 
Bi-2212^ and Bi-2223,^'^^ we consider here the minimal 
tight-binding model featuring a VHS and a d-wave BCS gap. 
We then perform least-square fits to determine the band 
parameters, the gap size Aq, and a scattering rate F. Fits 
to some representative spectra obtained on the as-grown 
sample are shown in Fig. 4. We observe that the model 
is able to follow the variety of measured spectral shapes 
with good accuracy. The amplitude of the BCS gap deduced 
from the fits is reported in Fig. 4b. The fits correctly yield 
a vanishing gap for single-peak spectra. We thus find that 
the single-peak spectra exhibit the same shape as the VHS 
in a two-dimensional tight-binding band. In general, the 
spectra measured in Bi-2201 are considerably broader than 
in Bi-2212 or Bi-2223. This hallmark can be captured by 
a phenomenological energy- independent scattering rate F. 
The typical values of F are indeed larger than in other Bi- 
based cuprates where F < 10 meV,^'^'^^ and we observe a 
tendency of F to follow the trend of the gap, as expected if 
both would result from the same pairing interaction. 



DISCUSSION 



FIG. 4. Theoretical interpretation of Bi-2201 tunneling spec- 
tra, (a) Comparison of representative dl/dV curves in the as- 
grown sample (circles) with fits based on a simple model (lines) . 
The model involves a two-dimensional tight-binding dispersion 
with a VHS, a d-wave BCS gap, and a constant scattering rate. The 
curves are offset vertically for clarity, (b) Squares: Amplitude of 
the BCS gap determined by the fits, as a function of — |£p|. When 
Aq = (single-peak spectra) the model corresponds to indepen- 
dent electrons subject to a constant scattering rate and the peak is 
the Van Hove singularity. Circles: Trend of the phenomenological 
scattering rate F as a function of — 



pseudogapped state to the un-pseudogapped normal state 
at fixed temperature (Fig. 3c). Since it is not possible to 
overdope Bi-2212 or Bi-2223 so as to reach a as low as 
in pure Bi-2201, ^'^^ such a transition has not been observed 
before by STM and can now be investigated because of the 
identification of the VHS. 



In Bi-2212 and Bi-2223, the bi- and tri-layer compounds 
of the same family, the spatial variations of the tunneling 
spectra have been ascribed to inhomogeneities of the doping 
level on a nanometer scale. In this scenario the single-peak 
structure we observe would correspond to extremely over- 
doped regions which are non-superconducting at 5 K. One 
thus expects that the proportion of such spectra increases in 
samples with higher doping levels. We have confirmed this 
expectation by studying overdoped samples with various 
doping levels. The STM observation of the VHS in a low-T^ 
cuprate superconductor provides a clear-cut experimental 
answer to a long-standing theoretical question. Our results 
confirm that the STM tunneling conductance is proportional 
to the full electron DOS, not only to the "superconducting 
DOS", and disqualify qualitative arguments used to claim 
the opposite. This must be taken seriously for the interpre- 
tation and modeling of STM tunneling spectra in materials 
with singularities in the electron dispersion. 
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The appearance of the VHS as a sharp peak in Bi-2201 
poses the question of its signature in the spectra of other 
Bi-based cuprates. It has been shown in Refs 5-7 that the 
couphng to the (n, n) spin resonance suppresses the VHS 
peak and gives rise to the strong dip in Bi-2212 and Bi-2223 
spectra. As suggested by ARPES^^ and our measurements, if 
a dip is present in pure Bi-2201, it must be very weak (see 
Fig. 4a) . Within a spin-fluctuation pairing scenario, we may 
tentatively ascribe the low of pure Bi-2201 and the ab- 
sence of dip in the DOS — hence the presence of the VHS — to 
a weak coupling to spin fluctuations. Interestingly in cation- 
doped Bi-2201, which exhibit higher critical temperatures, 
the VHS has not been seen.^^"^^ 

Our observation in the heavily overdoped sample (Fig. 3) 
confirms the presence of the pseudogap phase well into the 
overdoped region in Bi-2201.^ In a concurrent version of 
the phase diagram without a pseudogap in the overdoped 
region, the gapped spectra in Fig. 3b would have to be inter- 
preted as regions that are superconducting at T = 5 K > 
because of inhomogeneities. In this scenario, a mechanism 
different from the pseudogap would be needed to explain 
the anomalous large gap of ~ 30 mey 30 times larger than 
the expected value for a BCS superconductor with a of 
5 K. Furthermore the observation above of the coexis- 
tence of gapped and un-gapped regions strongly supports 
the existence of the pseudogap in the overdoped side and 
questions the universality of a cuprate phase diagram in 
which the pseudogap phase would cross the superconduct- 
ing dome. 



METHODS 

STS measurements. STS measurements were performed 
in a home-built STM placed in a low-temperature and ultra- 



high-vacuum environment similar to the one reported in 
Ref. 22. The crystals were cleaved at room temperature in a 
~ 10~^ mbar atmosphere, and immediately cooled down to 
5 K. The bias voltage was applied to the sample, while the 
ground electrode was a mechanically-sharpened iridium tip. 
The tunneling current was 0.25 nA and the regulation bias 
was 0.4 V Conductance maps were obtained by numerical 
differentiation of the I(y) curves. 

Modeling of the spectra. The dl/dV curves 
were modeled as the thermally broadened DOS, 
dl/dV oc j dco[-fXco - eV)]Nico), where / is the 
Fermi function and iV(co) is the BCS DOS given by 
^^M = Sfc(-lM)Im[a)-?fc+ir-|Afc|V(co+?fc+ir)]-\ 
The d-wave gap reads = Ao(cos/c^a - cos/c^a)/2, and 
for the dispersion we used the five-neighbor tight-binding ex- 
pression = 2ti(cos k^a + cos kya) + 4t2 cos k^acos kyQ + 
2t3(cos2A:^a + cos2kya) + 4t4(cos2A:^acos A:^a + 
cos/c^acos2/c_ya) + 4t5 cos2/c^acos2/c^a - /i with a 
being the lattice parameter. The k-mesh involves 
1024 X 1024 points. The model parameters were deter- 
mined by least-square fits to the STM spectra. Given the 
restricted energy range of the experimental data, the fits 
can not reliably determine the overall bandwidth controlled 
by ti. We therefore kept fixed to the ARPES value^^ 
ti = — 213 meV for all spectra. 
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